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We describe the growth of vesicles, due to the accretion of lipid molecules to their surface, in terms 
of linear irreversible thermodynamics. Our treatment differs from those previously put forward by 
consistently including the energy of the membrane in the thermodynamic description. We calculate 
the critical radius at which the spherical vesicle becomes unstable to a change of shape in terms 
of the parameters of the model. The analysis is carried out both for the case when the increase in 
volume is due to the absorption of water and when a solute is also absorbed through the walls of 
the vesicle. 
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I. INTRODUCTION 

Vesicles are small cell-like structures in which the mem- 
brane separating the contents of the vesicle from the en- 
vironment takes the form of a lipid bilayer. Part of their 
appeal comes from the fact that living cells are essentially 
very complex vesicles — with the membrane containing 
mixtures of different lipids and other components, a cy- 
toskeleton and complex surface structures [H . This has 
led to the use of vesicles as the basic component of models 
of protocells @, S On the other hand simple vesicles, 
without any additional structure, have many fascinat- 
ing properties when observed in the laboratory. Their 
self-assembly, their growth, their shape and the fact that 
they divide to produce daughter vesicles have many as- 
pects which are little-understood. The latter property 
of replication is especially interesting in the context of 
models of protocells. One can ask: how much of simple 
protocell dynamics can be explained using the statisti- 
cal thermodynamics of vesicles, without the introduction 
of more complex processes or of genetic material? This 
question will be the motivation for the present work. In 
particular we will be interested in describing the dynam- 
ics of vesicle growth and the instability which leads to 
the vesicle changing shape. 

To begin the construction of a model for vesicles, it 
is first necessary to review their dimensions. The typi- 
cal thickness of the bilayer is 7nm, while the radius of 
the vesicle itself can be anything from several times this 
value up to 100 /zm for so-called giant vesicles. There- 
fore the bilayer can be thought of as a thin membrane 
or shell enclosing the contents of the vesicle. In the case 
of biomolecules the standard picture of this membrane 
is the fluid mosaic model ^] where the proteins, enzymes 
and other such constituents are embedded in the lipid 
bilayer in which the lipid molecules can freely move as in 
a liquid. The proteins and enzymes may be able to move 
from the inner part of the bilayer to the outer part which 
is in contact with the environment. In the case of interest 



to us here, these biological aspects are absent and only 
the fluid nature of the bilayer remains. Therefore the pic- 
ture which emerges is of a two-dimensional surface, which 
is supplemented with a thin fluid layer on either side to 
describe the physical aspect of the bilayer Q . Turning 
this characterization into a quantitative description can 
be achieved in several different ways; the resulting mod- 
els go under names such as the spontaneous curvature 
model ff\ or the bilayer couple model [1, [13] . However 
these models have the common feature of a bending en- 
ergy of the membrane, which is given in terms of the 
curvature of the two-dimensional surface, together with 
some extra feature which, in some very basic way, ac- 
counts for the fact that the membrane has a thickness. 
In the spontaneous curvature model this is a extra fac- 
tor Co subtracted from twice the mean curvature of the 
surface, and in the bilayer couple model it is a constraint 
on the mean curvature. Here we will use the sponta- 
neous curvature model, for which the bending energy is 
E = {k/2) J dA {2H — Cof , where H is the local mean 
curvature, A is the surface area A, and k is the bending 
rigidity. 

The vast majority of studies of vesicles based on this 
form of the bending energy have been of a purely static 
nature; the energy E has been minimized subject to the 
constraints that the area of the surface. A, and the vol- 
ume of the vesicle, V , are kept constant. The first con- 
straint follows from the large elastic compression modu- 
lus; the energy scale associated with this is much greater 
than that associated with the curvature elasticity and 
this multiplies a term which fixes A Q. Similarly, the 
energy scale involving the osmotic pressure difference is 
so large, compared with the curvature energy, that the 
vesicle volume is effectively fixed. The shapes with the 
lowest bending energy are usually investigated at differ- 
ent values of a scaled spontaneous curvature and reduced 
volume. A "phase-diagram" of these minimal shapes in 
these two variables has quite a complicated structure, 
with "phases" of prolate ellipsoids, dumb-bell and pear 
shapes appearing, amongst others [6|. 
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Our intention here is to provide a means of linking 
these snapshots of the vesicle shape. To do this we need 
a dynamics which gives us a rule to move from one shape 
to another. We will be chiefly concerned with setting 
up the correct dynamical description of this system and 
so we will limit our attention to the growth of spherical 
vesicles and the transition from a spherical to an ellip- 
soidal shape. Our treatment differs from the previous 
studies of this problem [Tl|, [HI, in ways that we 
describe in detail in Section |IT1 However all these ap- 
proaches have the common features that they use macro- 
scopic variables and assume that the growth is sufhciently 
slow that it may be described within the formalism of lin- 
ear irreversible thermodynamics. 

The outline of the paper is as follows. In Section |TT] 
we introduce the formalism that will be used in the in- 
vestigation and compare our approach to those used pre- 
viously. In Section IIIII we carry out the analysis of the 
growth and loss of stability of the spherical vesicle in the 
simplest case of a purely aqueous environment and in 
Section IIVI we show how this generalizes when a solute 
is present. In Section |V] we summarize our results and 
discuss them in the light of previous work and the model 
assumptions. There is a mathematical appendix which 
gives the technical details relating to Sections Hill and HVl 



II. DYNAMICAL DESCRIPTION 

In the spontaneous curvature model the membrane is a 
two-dimensional surface, 5, which separates an inner re- 
gion, /, from the environment, E. The outer region could 
be a purely aqueous environment, or it could also contain 
a solute, with both the water and solute molecules being 
able to permeate through the membrane. 

The surface is a purely geometric construction; it con- 
tains no matter and simply has a bending energy asso- 
ciated with it. If its shape is known, then this bending 
energy only depends on the volume, V, it encloses: 



(Ci + C2 - Co)' dA, 



(1) 



where Ci and C2 are the principal radii of curvatures of 
the surface. In common with the other studies of this sys- 
tem [III, [13, [1^ [3 the deviations from equilibrium will 
be taken to be sufficiently small that the thermodynamic 
relation TdS — dE+PidV+'^^ fJ'adCa can be used. Here 
Pi is the pressure of the fluid inside the membrane, 
the chemical potential of chemical a which has concen- 
tration Cq. The only contribution from the membrane 
is a term dE = {dE/dV)dV , so changing the pressure 
inside the vesicle from Pi to (fj)cff — Pi + [dE/dV). 
Therefore, as long as we replace the internal pressure by 
this effective pressure, then we may ignore the membrane 
from a thermodynamic point of view, and simply treat it 
as a boundary which separates the inside of the vesicle 
from the environment. 



A. Purely aqueous environment 

The thermodynamic analysis of transport through a 
membrane which is simply a geometric transition region 
between two homogeneous regions was carried out by Ke- 
dem and Katchalsky fifty years ago [H, [l3|. Suppose, to 
begin with, that there is no solute present, so that only 
the flow of water through the membrane need be consid- 
ered. Then the usual assumption of the thermodynamics 
of irreversible processes, that the processes under con- 
sideration are sufficiently slow to giv e a linear relation 
between the fluxes and the forces [i3| , lead to the rela- 
tion [H, [3 



Jw — LpAP . 



(2) 



Here is the flux of water from the environment to the 
interior, Lp is the hydraulic conductivity of the mem- 
brane and AP is the difference between the exterior pres- 
sure, Pe, and that of the interior. However, as discussed 
above, to include the contribution coming from the cur- 
vature of the membrane we need to replace AP by the 
effective pressure difference given by 

(3) 

These results may now be brought together. The vesi- 
cle is assumed to increase its surface area due to extra 
lipid molecules being added to the surface. This in turn 
will change the pressure in the interior, and so change 
(AP)eff and give rise to a flux of water through the mem- 
brane. The rate of increase of the volume of the vesicle 
will by given by 



AP-(^ 
dV 



(4) 



using Eqs. ^ and ©. 

We now assume a growth law for the surface area, 
that is, the rate at which components are incorporated 
into the vesicle membrane. The simplest, and also the 
most plausible, is that this is proportional to the surface 
area 



— = XA 
dt 



A{t) = A{0)e 



At 



(5) 



The analysis below can be carried out with other growth 
laws. At a more fundamental level we would expect the 
correct form to emerge from the chemical reactions un- 
derlying this process. It is convenient to deflne a reduced 
volume by 



yl3/2 



(6) 



so that V — 1 for a sphere and u < 1 for all other shapes. 
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Then 



dv _ 3 6TT^/'^dV 
— _ --Xv + —^ — 

3, 6^1/% - 



AP 



dE 
dV 



(7) 



where we have used Eq. (jl]). 

If the spontaneous curvature, Cq, was absent in the 
definition of the bending energy Eq. Jl]), then the energy 
would be scale invariant: a typical length scale associated 
with the vesicle could be changed by an arbitrary scal- 
ing factor and E would remain unchanged. However, the 
inclusion of Co, which has dimensions of inverse length, 
introduces a scale into the problem. Suppose that R is 
the typical scale factor associated with the vesicle, then 
scaling the coordinates in Eq. ([T]) by this factor, the typi- 
cal length scale associated with the vesicle is unity. After 
rescaling we denote the principal radii of curvatures by 
ka = CaR, a = 1,2. Similarly, a dimensionless sponta- 
neous curvature, fco = CqR, may be introduced Q- If the 
vesicle is spherical, the typical length scale can be taken 
to be the radius, and Ci ~ C2 = l/R. Then 



E = 2ttk (2 - Coi?) 



dE 

w 



Cqk 

i?2 



(Coi? - 2) . (8) 



B. Including a Solute 



The formalism we have discussed in Section III Al can 
be generalized to include a solute. There will now be 
a flux of solute, J^, in addition to the flux of water 
Jw, and they will be linearly related to the thermody- 
namic driving forces which now includes the difference in 
the osmotic pressure of the solute across the membrane. 
Alls , as well as AP [l^, [l^ . The constants multiplying 
the forces in the linear relations are Onsager coefhcients, 
which will be symmetric in the usual way [17]. In fact 
we will use the linear relations involving a slightly differ- 
ent linear combination of variables, corresponding to the 
"second set of practical phenomenological equations" of 
Kedem and Katchalsky [16]: 



J„ = LpiAP-aAUs) , 
Js — c {1 — a) Jv + ujAHs 



(9) 
(10) 



where a is the reflection coefficient, c is the mean con- 
centration of the solute and ui is the solute permeability. 
The flux Jy is a linear combination of J^^ and , namely 
Jv = JwVw + JsV s where Vw and V s are the partial 
molar volumes of water and solute respectively. If we as- 
sume an ideal solute, then Alig = ksT Ac^^, where fc^ 
is Boltzmann's constant and Ac the difference in concen- 
trations across the membrane. This gives the more useful 
form 

J„ = Lp{AP - (jkBTAc) , (11) 
J, = c{l~a)Jy+ujkBTAc, (12) 



with Ac = Ce — Ci. 

In this case the total volume flow per unit area of the 
membrane is Jy, and so Eq. (jl]) is replaced by 



dV_ 
'dt 



Ajy 



LpA 



AP 



dE 
dV 



aksTAc 



(13) 



where, once again, we have replaced AP by (AP)off to 
account for the effect of the membrane curvature. Simi- 
larly, if N is the number of molecules of the solute in the 
interior, then 



dN 

— ^AJs = ^[c(l-cr) J„ 



- LjRTAc] 



(14) 



This last result may be written in a number of different 
ways using the relation N = CiV. 



C. Comparisons with Previous Work 

The analysis of the growth of vesicles presented in the 
next section will start from Eqs. (fT^]) and How- 
ever, we will end this Section by discussing how these 
two equations differ from those considered by previous 
workers investigating this problem. In Ref. |11| it was 
assumed that no solute was present, so the relevant dis- 
cussion is that presented in Section III Al The equation 
which was used was not Eq. ([4]), however, since the term 
involving the bending energy was introduced in a very dif- 
ferent way: the pressure difference in Eq. Q was simply 
set equal to (dE/dV), resulting in the pressure difference 
being completely absent from Eq. We believe our ap- 
proach to be the correct way of proceeding. The method 
of Ref. [ll[ is extended to include a solute in Ref. [13] ■ A 
further point of disagreement with our treatment is that 
the reflection coefficient, cr, is set equal to unity. How- 
ever this is only true if the membrane is impermeable to 
the solute, in which case to should equal zero too. The si- 
multaneous use of Eq. and Eq. (fH|) when a — 1 was 
already argued against in the original paper of Kedem 
and Katchalsky [l^] . 

The work reported in Refs. [13] and [13] has a differ- 
ent philosophy; there it is assumed that the instability 
by which the division process begins — the sphere be- 
comes unstable to an ellipsoid — is a Turing instabil- 
ity. Thus these authors introduce spatial effects, and 
the boundary of a two-dimensional vesicle is defined on 
a lattice. The "total pressure" involving the sum of the 
hydrostatic pressure difference, the osmotic pressure dif- 
ference and the term coming from the surface energy are 
all included. However, since this is a two-dimensional 
vesicle, the form of the bending energy is different, and 
it is not clear to us why the surface tension is included 
with such a large modulus. Thermodynamic relations 
similar to Eqs. p3|) and ^4]) are used, but apparently 
not to directly describe the change in vesicle shape. It 
is also not so clear to us what role the two metabolic 
centers that these authors introduce have in giving an 
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initial anisotropy to the vesicle and how crucial they are 
in initiating the symmetry-breaking instability. Clearly 
subsequent divisions will produce vesicles without these 
centers and therefore the mechanism will have to still 
work in their absence. 

One of the goals of the present work is to clarify the 
various assumptions made and systematize the method- 
ology that is used to study vesicle growth and division. 



III. THE FIRST BIFURCATION 

In the purely static analysis of the model defined by 
Eq. ([1]) , it is known that when the sphere becomes unsta- 
ble, the stable shape which replaces it is the ellipsoid . 
As an initial application of the formalism of Section [III 
we will investigate this instability from a dynamical view- 
point when no solute is present. 

The axisymmetric ellipsoid will be parametrized by ex- 
pressing the Cartesian coordinates as 

X = a sin 6 cos , 
y = a sin 6 sin <j) , 



X 



— c cos 6 , 



(15) 



where < < Stt, — 7r/2 < 9 < n /2 and where a and 
c are constants. For a sphere a = c = R, the radius. If 
the ellipsoid only differs in shape from the sphere very 
slightly, then a and c may be expressed as 



R{1 + aie) , c = i? (1 + cie) 



(16) 



where e is a small quantity and ai and ci are numbers 
which characterize the shape of the ellipsoid; if ai > ci 
it is oblate and if oi < ci it is prolate. 

Using standard results [lO] and p6)) , it is straightfor- 
ward to calculate the surface area and volume of the el- 
lipsoid for small e. The details are given in the Appendix 
where it is shown that 



l + -(2ai-Hci)e-f 0(e2) 



A = 4ttR^ 

V = ^^it!3 ^ (2ai +ci)e + (e^)] . (17) 

From these expressions we see immediately that the re- 
duced volume, V, defined by Eq. © is 1 + 0(6^), and 
so we have to go to next order in Eq. ([17]) to find the 
deviation of the reduced volume from the value 1 which 
it has when the vesicle is spherical. From Eqs. (|A4|1 and 
^51) we see that 



(18) 



with w < 1 for all cases except the sphere (ai ~ ci) 
as required. The bending energy is also straightforward 
to calculate, but much more tedious. From it we can 
determine dE/dV, as described in the Appendix. 



Substituting the expressions for v and dE/dV into 
Eq. © we find 
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\2 de 

(ai-ci) e—^Fo + Fie 
dT 



-0{e') , (19) 



where the Fj {j — 0, 1, 2) are functions of R, AP, Co and 
!]{= LpKC^ In2/A) and are given explicitly in Eq. (|A14[) . 
We have also introduced a scaled time r = t/Td, where 
Trf ~ In 2/ A is the time taken for the surface area to 
double in value. For Eq. (|19p to be consistent as e — > 
we require Fq = 0, which gives AP in terms of R, Cq and 
•q: 



AP 



C^R\n2 Cl 
2^ ^ "R 



2Co 
i?2 



(20) 



Since Eq. ([^O)) holds as e — > 0, it is true for the sphere 
and could have been obtained more directly from the 
bending energy of a spherical vesicle given in Eq. ([5]). 
Then, since for a sphere 



dV 
dt 



RdA _ XR 
'2~dt ^ ~2~ ' 



we see from Eq. Q that 



XR 



AP 



KCoRiCoR - 2) 
R^ 



(21) 



which agrees with Eq. ([^0)) . The quantity AP is the 
pressure difference required for the vesicle to remain a 
sphere while growing at a steady rate given by A. 

Setting Pq = in Eq. p^ . one sees that the two sides 
of the equation are not of the same order as e — > unless 
Pi = 0. From the explicit form for this function given 
in the Appendix, setting AP equal to the value given 
in Eq. ([^ gives the result displayed in Eq. (jAlSp . We 
see that, apart from making a particular choice for R in 
terms of Co and rj, we can only make Pi equal to zero 
by taking (2ai + ci) — 0. This simply amounts to a 
particular choice of ellipsoid shape. It is, in some sense, 
the most symmetrical choice and consists of changes in 
the direction of the two symmetric axes being half of that 
in the third direction (and having the opposite sign). 

Finally, setting both Po and Pi equal to zero, which 
implies the choice Eq. ([^0)1 and ci = — 2ai, gives the 
expression (IA17P for P2 and leads to 



de 
d^ 



( 1377 



977 



V8C2p2 



2C3P3 



13 In 2 
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(22) 



If the term on the right-hand side of this equation is posi- 
tive, then the sphere will be unstable to a transformation 
into an ellipsoid. This is the case if 277(13CoP — 36) > 
13CqP^ In 2, which can never be satisfied if CoP < 36/13. 
Similarly when viewed as an inequality which is cubic in 
CoP, one finds that it cannot be satisfied for any real 
positive CoP if 7y < 77min where 



3^ 2 In 2 
132 



17.94. 



(23) 
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For CqR > 36/13 the condition for the sphere to be un- 
stable may be written as 



1] > 



1 / l3Cgfl^ In 2 

2 \13CoR-36 



(24) 



The region of instabihty is shown in Fig. [TJ For values of 
Tj > rjmin a dynamical transition can eventually occur, the 
vesicle turning into an ellipsoid. We are only interested 
in the first transition which is encountered as R increases, 
which occurs for values of CqR greater than 36/13, but 
less than the value that corresponds to ry = r]^in- This 
latter condition gives the result CqR — 54/13. Therefore 
the critical radius, Rc, at which the transition occurs, 
lies in the range 36/13 < CqRc < 54/13. It is interesting 
to note the fundamental role that the phenomenological 
factor Co has in determining the critical radius. 



60 r 



40- 



20 - 





1 .• 

Ellipsoid 




\ . . ■ ■ 'm;n 


1 



10 



FIG. 1: The transition line defining the region of instabil- 
ity in the parameter plane (CqR, if) is depicted for a purely 
aqueous environment. There is no transition below the hor- 
izontal dashed line (which represents r^min = 17.94) and to 
the left of the vertical dashed-dotted line (which represents 
CqR = 36/13). In the remaining portion of the figure, the 
thick solid line represents the transition from the region on 
the left where the spherical configuration is dynamically fa- 
vored, to the region on the right, where the ellipsoid configura- 
tion is dynamically favored. The dotted line is an unphysical 
solution of the cubic equation (I26|l . 



So, in summary, we suppose that initially the vesicle 
is a sphere of radius _R(0). It then grows according to 
Eq. ([5]), that is, the radius increases according to R{t) = 
i?(0)e^*/2. The pressure difference between the interior 
of the vesicle and the exterior during this growth phase 
may be found from Eq. ([^0]) to be 



K 277 ' -R(O)' 



-Ai/2 



2Co 
i?(0)2 



-At 



(25) 

The growth phase continues until the vesicle has achieved 
a radius of Rc, which is the smallest real positive root of 
the cubic equation found by setting F2 equal to zero: 



liClRl In 2 



26r]CQRc + 7277 = . 



(26) 



As discussed above there are no such roots for 77 < r7min, 
given by Eq. and for rj > 77min, CqRc has to lie in the 
narrow range [2.77,4.15]. The critical radius is reached 
at a time 



t, 



ln{Rc/R{0)f 
In 2 



(27) 



At this time the spherical shape becomes unstable and 
the vesicle takes on an ellipsoidal shape. 



IV. THE BIFURCATION WITH A SOLUTE 
PRESENT 

In Section [n] we developed the formalism for the situ- 
ation where a solute was present, but for simplicity the 
analysis of Section IIIII assumed that the solute was ab- 
sent. In this Section we will repeat the analysis of Section 
IIIII with the solute included. 

The extra term which appears in Eq. (jl3p which 
changes the position of the instability is —aksTAc. To 
calculate it we need to determine Ac. This can be found 
from the other equation we introduced, Eq. (jl4p . however 
an integral over time has to be performed. To see this, 
we use Eq. to write Eq. (fH|) as 



dic,V) .dV 



(28) 



where c = Ac/ln(ce/ci) If we know A and V as 

functions of t, then we can in principle solve this differ- 
ential equation for Ci (t) and so find Ac(t) , which can then 
be substituted into Eq. (fT3|) . 

To illustrate these basic ideas, we will consider the spe- 
cial case where the membrane is impermeable to solute 
molecules, so that the number of solute molecules, N, 
is constant. However Ci{t) does change with time, due 
to the fact that the volume, V , increases with time, and 
this has a non-trivial effect on the instability analysis. A 
membrane impermeable to solute molecules is defined by 
a = 1 and a; = 0. It follows directly from Eq. (pS)) that 
CiV is a constant, and so 



C^it) 



V{t) 



(29) 



Then the term —ksTAc which appears in Eq. gives 
rise to an additional term on the right-hand side of 
Eq. ([7]) which equals 



6^1/2 



Lj^keT {ceV{t) - cMVm 



Using the expressions p7|) . this becomes 
iLpksT 



(30) 



1 



i2ai + ci) e + O U') 



i?4 l^* 3 
X {ceR^ [I + (2ai + ci) e + O (e^)] - c,(0)i?3(0)} , 

(31) 
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to first order in e, since initially the vesicle is a sphere. 
To see how these changes affect Eq. (|19p let us write it 

as 

— ^(ai-cif e^ = Go + Gie+G2e2 + 0(e3) , (32) 
15 In 2 dr 

where the Gj {j = 0, 1) are related to the Fj as follows: 

3r/$ [3r - AirR^Ce] 



Go — Fq 



Gi = 



inC^R^ In 2 
r/$ [3r - ttR^c^] 
ttC^R^ In 2 



(33) 
(34) 



where T = Ci{0)V{0) and $ = akBT. The analogous 
result for G2 is given by Eq. (jAlSp in the Appendix. Fol- 
lowing the same line of argument as in Section IIIII we 
require Go and Gi to be zero in order for the stability 
analysis to be applicable. The first condition gives an ex- 
pression for the pressure difference across the membrane: 



AP 



G^i? In 2 



+ 



2?7 



G2 

R 



2Co 

i?2 



{R^C- R^{0)c,{0)) 



(35) 



which shows the additional terms that are added to 
Eq. (POI) when a solute is present. The second condi- 
tion, Gi = again implies that (2ai +ci) = and so the 
addition of the solute does not change the shape of the 
ellipsoid for which the stability analysis applies. 

If we now use the conditions found by implementing 
Go = and Gi = in Eq. (|A18|) and substituting this 
into Eq. ([32]) we find 



de 
dr 



1377 
45r$77 



327rG(*i?4 



977 



2G3i?3 



13 In 2 
16 



(36) 



which gives the required modification of Eq. (I22|) . A sim- 
ilar analysis to that given in Section IIIII shows that if 
47rGoi?(13Goi? - 36) < 45r$, then the sphere is always 
stable, no matter what the value of 77. This gives the 
minimum value of the radius, R, which can lead to an 
instability (corresponding to an infinite value for rj) to 
be 



CoRn 



18 
13 



1 



1 



65r$ 

1447r 



If i? > i?min, the spherical shape is unstable if 
V > 



26TrC^R^ In 2 



47rGoi? (13Goi? - 36) - 45r$ 



(37) 



(38) 



Figure [5] depicts the domain of instability in the pa- 
rameter plane (Goi?, ?y). Different curves refer to differ- 
ent choices of the quantity r<I>, and allow the qualitative 



inspection of the modification induced by the presence of 
a solute. The condition for a double root of the quartic 
equation for Goi? gives, as in Section IlIII the maximum 
value of R which can lead to an instability: 



27 

Goi?niax — Yg 



65r$ 

1627r 



(39) 



with the corresponding value of rj being rjmin- This can 
be determined as a function of r<I> by substituting the 
expression (j39p for Goi?max back into the quartic equa- 
tion. The resulting function is shown as an inset in Fig. [51 
This increases as r$ increases (larger solute concentra- 
tion and/or higher temperature). However, the ratio of 
-Rmax/^min remains essentially unchanged at about 1.5 
for all values of r$. This means that the range of val- 
ues of the radius at which the spherical shape becomes 
unstable remains quite small. 




FIG. 2: (Color online) The region of instability is depicted for 
the case when a solute is included in a similar manner to that 
shown in Fig.[T]when a solute is not present. Different curves 
refer to distinct values of the product r"I> {— 0, 5, 10, 20). The 
dotted line which passes through the minima of these curves 
shows the maximum value of CqR for which a transition can 
occur, and is given by Eq. (|39|l . The corresponding value of 
77, T/min, is plotted in the inset as a function of the quantity 

r<i>. 



V. CONCLUSION 

Despite the great interest in the growth, change in 
shape and division of vesicles, very little is known about 
the nature of the processes that govern them. Even fun- 
damental questions about the typical vesicle radii at the 
various stages or the hydrostatic or osmotic pressure dif- 
ferences between the exterior and interior, are still largely 
open. The main obstacle to achieving a greater under- 
standing is the difhculty in carrying out experiments. 
Even qualitatively, a consistent picture is hard to achieve, 
and so a theoretical description which would help with 
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the interpretation of experimental results would be very 
welcome. In this paper we have proposed such a de- 
scription, taking extra care to correctly incorporate the 
energy associated with the curvature of the membrane in 
the thermodynamic description. We have concentrated 
on establishing the formalism and demonstrating it on 
the initial stages of the growth and on the first bifurca- 
tion to an ellipsoid. The subsequent time evolution of 
the vesicle, leading to division, can now be investigated, 
but this can only be carried out numerically and we leave 
it for the future. 

The vesicle growth is ultimately caused by the incor- 
poration of lipids from the environment into the vesicle 
wall. This will increase the area, which in turn will lead 
to a change in the internal hydrostatic pressure and so 
initiate a flow of fluid through the membrane. This will 
happen extremely slowly — at a rate governed by the 
parameter A, and so in practice the vesicle will never 
become flaccid and will maintain a spherical shape in 
the initial stages of the evolution. It is the quasi-static 
nature of the expansion that allows us to use the for- 
malism of linear irreversible thermodynamics p^j . This 
predicts, among other things, that the radius of the vesi- 
cle will grow according to R{t) = i?(0)e^*/2 = 2*/(2Td) 
using A = In2/Tcj. However this exponential growth gets 
cut off at a radius, Rc, given by the smallest real solu- 
tion of Eq. at a critical time given by Eq. ^7\ . This 
critical radius falls into a remarkably narrow range. The 
corresponding value of the hydrostatic pressure can then 
be determined. 

It is interesting to compare these results with those 
found by the previous studies discussed in Section Hi CI 
In Refs. [Ill, [I4I , the term involving the bending energy 
on the right of Eq. Q was effectively omitted and there- 
fore the expression pp)) did not contain the final two 
terms, which originate from this term (see Eq. ([8])). It 
could be argued that it is a legitimate approximation, in 
the sense that for the parameters of interest, this term 
is negligible compared to the one that is retained. Since 
we have included both contributions, we are able to test 
this by calculating the ratio between the first term on the 
right-hand side of Eq. (PO)) and the sum of the second and 
third terms. We find that this ratio always lies between 
and 1. That is, the terms omitted in Refs. [H!, [13 are al- 
ways more important than the one included. This ratio is 
actually zero at one extreme of the range of allowed val- 
ues of R {CqR = 36/13, which corresponds to an infinite 
value of the parameter 77), and increases monotonically 
to 9/14 when CqR takes on its greatest allowed value 
(54/13) and 77 takes on its smallest allowed value (17.94). 
It should also be remarked that in Rcf. 11] the value of 
r] used is 1.85, which is the value at which vesicle split- 
ting gives rise to two equal vesicles, which are equal in 
size to the initial vesicle. Given the inconsistencies in the 
formulation we have just alluded to, it is doubtful that 
this value is correct. Further work is required to enable a 
comparison between the work reported in Refs. [13] and 
[T3 ] and our approach. As explained in Section Hi C[ the 



method used is very different: the vesicle is a closed curve 
in two dimensions, metabolic centers are present which 
induce the symmetry breaking and the results obtained 
are purely numerical. 

These results, and others involving the nature of the 
initial bifurcation, provide an initial set of predictions 
which can be compared with future experiments. Of 
course, the model as it stands is quite simple and many 
features could be made more realistic. For instance, Cq 
is almost certainly not a constant, and becomes non-zero 
during the course of the vesicle growth. Nevertheless, 
taking a non-zero Co induces a transition [2l| . Another 
aspect that we have not included in the present treat- 
ment is a discussion of the effect of thermal fluctuations 
on the shape of the vesicle. At first sight it is not obvious 
if thermal fluctuations will have any appreciable effect. 
An order of magnitude calculation suggests that un- 
der most circumstances fluctuations will not need to be 
taken into account, however near an instability there may 
be a significant effect. Within the static picture, which 
has been the subject of by far the majority of papers to 
date, the procedure for inve stig ating the effect of ther- 
mal fluctuations is clear [l^, [231 f24l . |25| . |26| . The shape 
is written as a sum of a stationary term plus a small 
fluctuation and the energy expanded to quadratic order 
in the fluctuation. Putting this into a Boltzmann factor, 
the Gaussian integrals may be performed. However there 
are a number of technical issues, such as the inclusion of 
constraints. There is also the question of the size of non- 
Gaussian fluctuations. Nevertheless, it is clear that ther- 
mal fluctuations can shift the position of an instability, 
or even change its nature completely. The analogous cal- 
culation carried out within the present framework would 
be even more complex, however it should definitely be 
addressed in future work. 

Despite such shortcomings, and there are undoubtedly 
others, we believe that the model is sufficiently detailed 
to provide a reasonably good description of vesicle growth 
and hope that it will serve to clarify a number of aspects 
in this fascinating, and neglected, field. 
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APPENDIX A: SURFACE AREA, VOLUME AND 
BENDING ENERGY OF AN ELLIPSOID 

In this Appendix we will collect together the results 
which are required to calculate the reduced volume and 
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the bending energy of an ellipsoid in Section imi together 
with the generalization in Section HVl 

The ellipsoids we consider here are axisymmetric. The 
prolate version is formed by rotating an ellipse with semi- 
minor axis a and semi-major axis c (i.e. c > a) about the 
major axis. It sm'face area is given by [20l | 



A = 27ra^ ■ 



2'Kac 



For an oblate ellipsoid, c < a, and [20|, 



In 



a + \/ cp- — (? 



In both cases their volume is given by [2C 



Y = —na c . 



(Al) 



(A2) 



(A3) 



Substituting the parametrizations (|16|) into the expres- 
sions (IA1[) - (IA3|) one finds Eq. pT|) . As mentioned in the 
text, this immediately implies that the reduced volume 
is 1 to this order, and so we have to go to next order in 
e to see some deviation from the result for a sphere. At 
next order 



A 



1 



1 + -{2ai + ci)e 



+ — {6al + cl + Saici) + O (e^) 

±0 



(A4) 



and 



V 



^nR^ [1 -I- (2ai +ci)e 



(A5) 



which together give the expression (|T8|) for v to second 
order in e. 

The other quantity we have to evaluate is the bending 
energy, E, given by Eq. H]). This involves evaluating the 
two integrals 



Ji 



HdA, Jo 



dA. 



(A6) 



where H = (Ci + C2)/2 is the mean curvature. For an 
axisymmetric ellipsoid this is given by [20| 



H 



c [3a2 + + (a^ - c^) cos 20] 



cos^ 



c2sin2 ef^^ 



(A7) 



Evaluating the integrals in Eq. (jA6|) using the result (|A7I 
yields 



Ji = 2Trc + 
Ji = 2ttc + 



c + \J (? — a? 



27ra2 



Va2 



(c > a) 
(c < a) , 



(A8) 



and 



J2 



TTC 

a 



7a 2a^ . , 



J2 = 





'tanh ^ 






a 




Va2 


-C2 



7a 2a3 

3c2 ^ 3c4 



(c < a) . 
(A9) 



Substituting the parametrizations (|16p into the expres- 
sions (|A8|) and l|A9[) . and using Eq. (|A4p . one obtains 

= 27rK j(Coi?-2)^-l-^(2ai-hci)Coi?(Coi?-2)e 



1 

15 



[c? (32 - 4Coi? + C^i?^) + 8aici (-8 + Coi? 



ClR^) + a? (32 - 4Coi? eCo^i?^)] ^2 _^ ^^3j | 

(AlO) 



The aim is to calculate dE j dV which we achieve by using 

(All) 

From Eq. (jAlOp one finds 



dE _ dE dR 
dV ^ 'dR'dV 



— = 47r^Co|(Coi?-2) + -(2ai-f ci)(Coi?-l)e 

+ [cl [CoR - 2) + 4aici (1 + 2CoR) 
15 

+ al{6CoR-2)]e^ + 0{e'^)} 
and from Eq. (|A5|) 



(A12) 



= 47ri?^[l + (2ai + ci)e 
+ (a2 + 2aici)e2-^0(e3)] , (A13) 



which gives an expression for dE/dV. Substituting this 
into Eq. ([7]) , and making use of the expression (|18p for v 
we obtain Eq. (fT9|) with the i^j {j = 0, 1, 2) given by 



^^0 - 3 



2 - 4Cor/ -I- 2Q?7i? - 2APryi?2 
2C(^i?3 In 2 



i^i = (2ai-Hci)?? 
1 



i^2 - 



15C^R^ In 2 



( 6Co - 2C^R + APR^ \ 
\^ C4i?3 1n2 j ' 

[2aiCi(-82Co77-(- 16C^r/i? 



- 9APriR^ + 6ClR^ In 2) - c? (98Co?7 - 29C^r]R 
+ QAPrjR^ + 6C(Ji?3 In 2) - al (278Co?y 

74C^77i? + 2lAP77i?2-^6C^i?Mn2)j , (A14) 

where AP = AP/k. 
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As explained in the main text, setting Fq = gives 
Eq. ([^0]) . Using this expression for AP then gives 



Fi = (2ai +ci) 



^C^R^ ln2 + St] ~ 2riCQR 



ClR^ In 2 



(A15) 



Apart from exceptional values of R, this vanishes only 
when 2ai + ci = 0. 

If we now set ci = — 2ai we find from Eq. (jA14[) that 

- 1^ 5C4i?3in2 j • 

(A16) 

Finally, setting AP to the value given in Eq. ([20|) one 
finds 



^ „ 2 /^26ryCoP- 13Cf^p3ij^2-7277 
-^2 — oa^ I 



lOC^pa In 2 



which leads to Eq. ([22 



In Section IIV| the calculation is carried out in the 
present of a solute. The equation which describes the 
instability is now Eq. p2p with Go and Gi given by 
Eqs. ([33)1 and (p4|) respectively and G2 given by 



G2 — F2 — 



77$ 



107rGi^p4 



[4aici (2ir - 37rP^Ce) 



+ cl (33r - 47rp3c^) + 2a? (54r - luR^Ce)] , 



(AI8) 



(A17) where once again T = Ci{Q)V{Q) and $ = aksT. 
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